Baryon number conservation and enforced electric charge neutrality for bulk viscosity 

in quark matter 
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General constraints on fluid velocity divergences for particles in quark matter are derived from 
baryon number conservation and enforced electric charge neutrality. A new oscillation pattern 
in three-flavor normal quark matter satisfying these conditions is found and its bulk viscosity is 
calculated. The result may have astrophysical implication for maximum rotation frequencies of 
compact stars. 
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In the core of compact stars where the density could 
reach 5-10 times the normal nuclear matter density, the 
constituents of nucleons and hadrons can be squeezed 
out to form quark matter [l|, Q • Quark matter in normal 
state could exist in various forms, such as strangelets 
0, 0, 0], mixed phases @, 0, S, S| and etc.. Search- 
ing for quark matter in stellar objects is very challeng- 
ing. The cooling behavior dominated by neutrino emis- 
sions may be able to distinguish nuclear and quark mat- 
ter pH EES El El EI EE EE ■ But it is still sub- 
jected to many experimental and theoretical uncertain- 
ties firl . Esl . [l9| . Other transport properties such as 
shear and bulk viscosities are also of great interest to 
this end. The shear viscosity damps the differential ro- 
taion to make a uniform rigid body rotation of stars. The 
bulk viscosity is crucial for the damping of pulsations in 
compact stars. Such pulsations could take place during 
the formation of stars or by external perturbations. They 
could also be driven by instabilities of gravitational wave 
radiation such as the r-mode instability which arises in 
the absence of damping effects. There has been a lot of 
literature about calculations of viscosities in nuclear mat- 
ter in various situations |2(| EH, E2, 0, EI HE EE . The 
shear viscosity for color-flavor-locking phase has been cal- 
culated [13] ■ The bulk viscosities for normal and color 
superconducting quark matter have been studied exten- 
sively [II, El EE OS El El- In this letter we put empha- 
sis on general constraints of baryon number conservation 
and charge neutrality on the bulk viscosity. Following 
these constraints, we will give new solutions to the bulk 
viscosity in normal quark matter different from previous 

results (U EI EE EE Ei Ei El . 

The bulk viscosity is associated with the damping 
of the baryon density oscillation denoted by Sub = 
8riBo eZUJt with the amplitude Shbo and the frequency w in 
the range 10 3 — 10 4 s _1 . We assume that the amplitude 
of the oscillation is small and can be treated as perturba- 
tion to its equilibrium value. The perturbation in baryon 
density drives quark matter out of chemical equilibrium 
via following processes 



u + d «-» u + s, (d <-> s transition), 
u + e «-> d + v, (Urea I) , 
u + e «-» s + v, (Urea II). 



Here we consider the normal quark matter with three fla- 
vors u, d and s. Light quarks u and d are almost massless, 
while s quarks have large mass. The relaxation to chem- 
ical equilibrium is related to deviations of chemical po- 
tentials from their equilibrium values for these processes, 

6 Hi = Hs - Hd = $Hs ~ SfXd, 

6^2 = Hd - Mu ~ He = <W - ^« ~ <fy«e, 

8 Hz = Hs - Hu - Me = Shs - o~Hu - 8 He- (2) 

The reaction rates for processes in HJ can then be written 
as \k8nk (k = 1,2,3) for small Snk <SC T, where Xk are 
coefficients. In our calculations, Ai and A3 are taken 
from Ref. (33l . E3| except that we have included the mass 
effect of s quarks [3(3] and phase space reduction from 
Fermi liquid behavior [T3, E3, EE- The value of A2 is 
taken from Ref. 

With massive s quarks a purely three-flavor system 
cannot be electrically neutral, so there must be electrons 
in the system. In reactions (JTJ) , the electron number and 
flavors are not conserved since they can be created or de- 
structed, however the baryon number and electric charge 
are conserved. The continuity equations for quark flavor, 
baryon and electron number read 



dm 

— + UiV ■ Vi 
dt 



Ji, 



u, d, s, B, e) 



(3) 



where the substantial or material derivative is defined by 



— d 



T£ + Vj • V. The sources Ji are linear combinations 



of reaction rates \k8nk when Snk are small. It is obvious 
that Jb = required by baryon number conservation. 
The following relations for baryon number conservation 
and charge neutrality are widely used in literature, 



riB 



i—u,d,s 



E ■ 

i—u,d,s 



Jb = \ E j * = < Je= E ®* J i- ( 4 ) 



i—u,d,s 



i—u,d,s 



(i) 



The charge neutrality condition in the above is a natural 
constraint since the accumulation of net charges would 
make the Coulomb energy density blow up and be less 
favorable in energy. With Eq. and continuity equa- 
tions ([3]), one obtains following constraints on number 



2 



densities and velocity divergences, 

n B \7-v B = ^2 g^V-v,, 

i— u,d,s 

n e V • V e = QiUiV ■ V,; 

i— u,d,s 



(5) 



where the first line is the result of baryon number conser- 
vation and the second one that of charge neutrality. Eq. 
([5]) was not presented and applied explicitly in previous 
literature. A simple solution to the above equations can 
be found at first glance [H, \M, E3, El, S] , 



V • Vj = V-v B , (i = u, d, s,e), 



(6) 



with number densities satisfying Eq. This solution 
corresponds to the homogeneous or one-component fluid 
with a single fluid velocity field. Inserting Eq. © back 
into continuity equations J3]), one finds 



n B - 



dt 



drii 
~dt 



Y dTlB - T 
-A-i — TT — J ii 

dt 



(7) 



where Xi = rn/ns are partial fractions for particles i in 
terms of the baryon number density. We have used the 
continuity equation for the baryon number in Eq. JH}, 
Obviously Eq. J7]) respects baryon number conservation 
and charge neutrality. A property of the solution is 
that it will give a vanishing bulk viscosity in the case 
that all quarks are massless, where the system is always 
in chemical equilibrium. An example of this property in 
two-flavor normal quark matter can be found in Ref. [33j]. 

In this paper we rigorously apply the constraints J5]) 
and find a new solution to Eq. © different from Eq. . 
Our starting point is that the role of strange quarks is 
special since they are much heavier and may respond to 
the density oscillation more reluctantly than other parti- 
cles. As an extreme case, we assume V • v s = 0. We also 
assume V • v n = V • v^. Then one obtains from Eq. (O 
the following relations, 



V • v e 

V • V u .d 



n B (2n t , 



n dj 



n e (n u 
n u + n d 



4- rid) 

■V-VB 



V-VB, 



(8) 



Generally the velocity divergence for a particle in a fluid 
depends on its mass. Rigorously the fluid with more 
than one particle species with different masses should 
be treated as a multi-component fluid. The above solu- 
tion to Eq. ([U is reasonable in the case that the masses 
of strange quarks are of the same order as the chemi- 
cal potentials and much larger than light quark masses. 



If the strange quark masses are small, one can inves- 
tigate many other solutions which are close to the so- 
lution JSJ, for example, one can assume velocity diver- 
gences of particles deviate from that of baryons by a small 
amount, V • v» = V • + £j with (i = u, d,s,e). Fol- 
lowing Eq. |(5]), these e, satisfy J2i=uds \ n i £ i = and 
a s Qi n i £ i = n e£e- Any small values of e, under 
these constraints denote a slighly different solution from 
Eq. ([6]). We will not consider these solutions in this pa- 
per and focus on the solution (JSj) - With these relations 
in Eq. ([8]) for velocity divergences the bulk viscosity for 
a system of quarks and electrons can be derived. The 
deviation of the pressure from its thermodynamic value 
is related to the bulk viscosity [35| . 



SV 



-CV-VB. 



(9) 



The above also provides a definition for the bulk vis- 
cosity. The variation SV can be expressed in terms 
of density variations Srii = SniQe tuJt for some quarks 
or electrons, which are linearly independent after using 
Eq. ((!]). Here Sriio are complex amplitudes and can be 
solved by applying continuity equations. The number 
of continuity equations applied is equal to that of inde- 
pendent densities in order to close the system of equa- 
tions. Normally the baryon density n B is set to be one 
independent variable. The r.h.s of Eq. ([9]) is actually 
(,(n B )~ 1 dn B I 'dt. So a complex £ can be obtained from 
Eq. ©, C = n B SV/(dSn B /dt), whose real part gives the 
bulk viscosity [35l |. 

We will first consider in reactions two simplest 
cases, the d <-» s transition (by turning off Urea I and 
Urea II) and Urea II (by turning off d <-> s transition and 
Urea I) , separately. Finally we will address a more realis- 
tic case, the three coupled processes together. The Urea 
processes with light quarks (Urea I) was already studied 
in Ref. The calculations for the d «-> s transition 

without charge neutrality can be found in Ref. [23]. We 
list all quantities needed in evaluating bulk viscosities in 
three cases in Table [H The second column is for indepen- 
dent variables. The third one is for continuity equations 
used in this paper to solve the variations of densities. 
For example, we use continuity equations for s quarks 
and electrons for the three coupled processes. With Eq. 
([8]) it makes no difference to use other independent conti- 
nuity equations. The fourth column lists other variables 
expressed in terms of independent ones. In the d <-* s 
transition and Urea II, there are additional constraints 
on electrons and d quarks respectively, since they do not 
participate in the reactions. With J ejS = in Eq. J3]) 
and (]8]), one can solve Sn e and Sn d in terms of Sn B in 
the last column. 



I 

For the d <-> s transition, following Eq. |(9]) and with the second row in Table HI one obtains the bulk viscosity, 



( n B 5V \ Aingjyi 

\d8n B /dt) w2 + A 2^y 



Ci = Re 



Table I: Quantities for bulk viscosities in three cases with C = 3n t j/(n„ + n<j). 





independent variables 


continuity eqs. 


other variables 


special constraints 


d <-> s 


8n B , Sn s 


s 


5n u — (3 — C) Sub 
Srid = CSriB — Sn s 


Je = 

Jn e = (2 - C)<5n s 


Urea II 


Sub, 8n s 


s 


8n u = (3 — C)5tib — <57i s 
<5n e = (2 — C)5ub — <5n s 


Jd = 
5rid = CStlb 


coupled 


Sub, 8n a , 8n e 


s, e 


<5n u = Sub + 8n e 

5nd = 25ns — <5n e — <5n s 





Figure 1: The bulk viscosity £i for d <-> s transition. The 
masses and chemical potentials (in MeV) are set to m u = 
m d = 0, m s = 100, (j, s = ^ d = 500, = 495, ^ e = 5. 
Note that these values satisfy charge neutrality and chemical 
equilibrium conditions. 
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where 



dm _ 



dn B 

dm 

dn s 
dV 
dn s 



= -C 



dfi ti 



dn d 

dnd dn s 
dfi s 
an s 



n d 



djM 

dn d ' 



with C = 3rid/(n u + nd). One can verify that Ci is def- 
initely positive since < and ^ a < using the 
equation of state for degenerate Fermi gas for d and s 
quarks. The numerical results for Ci are shown in Fig. [TJ 
As shown in the figure that the bulk viscosity increases 
with decreasing frequency until it satuarates below a crit- 



ical frequency u>i c ~ Ai 



<9n. s 



One also sees that for low 



frequencies ui <C u>i c , £i is inversely proportional to the 
transport coefficient Ai. This means the faster the reac- 
tion proceeds the smaller the bulk viscosity is. 

Similarly, with the third row in Tab. HI the bulk vis- 
cosity for Urea II is, 



Figure 2: The bulk viscosity £3 for the Urea II. The parame- 
ters are the same as in Fig. [l] 
34 




where 



dfi3 
dn B 
dfi 3 
dn s 
dV 
dn s 



-(3-C)^p-(2-C) 



d[i s 
dn H 



dn u 
dfi u dn e 
dn v dn e 



dfi e 
dn R 



dfi u 
dn u 



d[is 
'<9n s 



dfi e 
dn e 



One sees §r~ 



< due to C < 2, and 



< 



dn s 3/i s 3/i 

with jji e <c jjL u and m s 3> m u - The numerical results 
are shown in Fig. [2j The behavior of £3 is similar to £1 . 



Sua 



< U>lc 



For critical frequencies one sees uj^c ~ A3 

because A3 < Ai. The saturation value of C3 is larger 
than that of £i for the same reason. At high frequencies 
u> ^> lusc one observes (3 <C Ci • 

With the fourth row in Table HI the bulk viscosity for 
three coupled processes can be expressed as 



c 



riB 



n t Srls 
D s lm- — 

u> \ ons 



n B 



D, [in— i = ~ ( D„N s +D e N e ), 

(11) 

where F and N e>s are quadratic and linear functions of w 2 
respectively. They also depend on coefficients A&. The 
explicit form of these functions will be presented else- 
where [36|| . The coefficients D e s are given by 
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Figure 3: The bulk viscosity £ for the coupled processes at 
T = 0.1 MeV. The dashed and dash-dotted curves are the 
results from the solution Eq. with two sets of parameters. 
The solid curve is the result from the solution Eq. ([8]). The 
parameters of the solid and dashed curves are the same as in 
Fig. [T] and [21 while those of dash-dotted curve are different 
and still satisfy neutrality and chemical equilibrium condi- 
tions, m s = 300, /i s = [id = 500, fi u = 456, /i e = 44, all in 
unit MeV. 




icr 3 i(r 4 i(r 3 icr 2 io-' io° io 1 io 2 



T[S] 

D s = n s - n d - — . 

on s On d 

The numerical result for the bulk viscosity £ is shown 
in Fig. [3] by solid curve. Also shown are dashed and 
dash-dotted curves from the method with Eq. J6]). With 
same parameters, the result from Eq. (J8j (solid curve) 
is about one order of magnitude larger than that from 
Eq. iJH) (dashed curve). When m s is set to 300 MeV, the 
result from Eq. ((H (dash-dotted curve) is comparable to 
the solid curve at high frequencies but still differs from 



it substantially at low ones. Generally there can be up 
to 3 plateaus in the log-log plot of the bulk viscosity as 
function of frequency depending on parameters. In Fig. 
03 the solid and dash-dotted curves have two plateaus, 
while the dashed curve has only one. It can be verified 
that C is dominated by the s quark part in Eq. (fTTjl . 
Since D s is negative, the phase of 5n s should be delayed 
relative to Sub manifested by the positivity of the bulk 
viscosity. The another calculation of the bulk viscosity 
for the coupled processes can be found in Ref. 

In summary, we have derived general constraints on 
fluid velocity divergences of quarks and electrons in nor- 
mal quark matter from baryon number conservation and 
enforced electric charge neutrality. Under these con- 
straints we find a new solution to velocity divergences 
leading to different bulk viscosities. As an extreme case, 
this new solution could be realistic in some circumstances 
where strange quarks respond to the oscillation in a very 
different way from light quarks. Other solutions are also 
allowed by these constraints. The new result for the 
coupled processes may have astrophysical implication for 
larger maximum rotation frequencies of compact stars. 
Similar constraints on fluid velocity divergences can also 
be obtained for nuclear matter. 
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